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An Effective Field Theory Model for One-Dimensional CH Chains: Effects at Finite 
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In this work we use an effective field theory model to investigate doped (CH)x chains under the 
influence of an external constant Zeeman magnetic field Bo, at zero and finite temperature, in the 
mean-field approximation and beyond. We consider both homogeneous and inhomogeneous A{x) 
condensates and calculate the Pauli magnetization and the magnetic susceptibility of these chains at 
various situations for the temperature and chemical potentials. We also briefly discuss the possibility 
of using these materials as partially polarized ID organic conductors. 
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I. INTRODUCTION 



In the last few years graphene has attracted the attention of both theoretical and experimental communities due 
to the interesting features this two-dimensional system presents. The electrons of this material are Dirac-like, with 
^ linear dispersion relations near the K points. Monolayer graphene can be cut in stripes to construct one-dimensional 

d (ID) graphene wires, which have been studied in various recent works [l], [l]- In Ref. Q, an effective field theory has 

. been used to investigate the low-energy quantum electrodynamics of Dirac electrons in an undoped graphene wire. 
I ' Another organic material which presents similarities with a graphene wire is a CH chain, which for more than 3 
' O , decades have also created an admiration due to the unexpected and fascinating discover that doped trans-polyacetylene 
C ' (TPA), a ID polymer, behaves as a metal, exhibiting electrical conductivity of some metals, like copper Since 
then, several promising properties, such as electronic, optical and magnetic, have been shown by conjugated polymers, 
which give these materials the possibility of application in the semiconductor nanotechnology [1, @] ■ 

Polyacetylene is a linear chain of CH groups which can have two forms, trans and cis. The trans form, is more 
stable and, in the absence of doping, has a doubly-degenerate ground state. The trans- (CH)x has one tt— electron per 
site, and would be a metal in the tight binding approximation. However, spontaneous Peierls dimerization (i.e., the 
interaction of the electrons with the lattice) turns TPA into an insulator. After the critical doping, the conductivity of 
TPA is greatly increased, reaching metal-like properties Q- Experiments show that the observed non-metal to metal 



, first-order transition in polyacetylene happens when the dopant concentration y, defined as the number of doped 



electrons per 100 carbon atoms, is increased up to a critical value — 6 [8|. Recent calculations show that thermal 
[ effects have a weak effect on yc [9] . 
■ The electron-phonon interactions in CH is successful described by the Su-Schrieffer-Heeger (SSH) Hamiltonian lo!], 
T— I ' whose continuum version is known as the Takayama-Lin-Liu-Maki (TLM) model [llj . The TLM model is a relativistic 
. field theory with two-flavor Dirac fermions. We shall use the Gross-Neveu (GN) model [l^ in 1+1 dimensions since 
. ^ ' its Lagrangian is equivalent to that of the TLM model Ell- The GN model is a well-known ID theory with four- 
, fermion interactions. The choice of employing an effective quantum field theory model (QFTM) to the investigation 
H ' we do in the present work is basically motivated by the analytical and simple solutions the GN model provides, as 
. y .' we will see below. Besides, one of the most important features of the GN model which enables its use as an effective 
theory to describe materials that undergo transitions from non-metallic to metallic phases, as the polymers we are 
considering, is that it is asymptotically free, just like QCD. This means that the fermions become non-interacting at 
high enough densities and are not able to support their condensate that they formed to break chiral symmetry at 
zero or very low densities. Therefore, since the theory displays a phase transition at a certain fermion density or at 
the corresponding value of chemical potential, the use of the GN model is very appropriate in this context. In other 
words, or in the field theory language, the discrete chiral symmetry is spontaneously broken in the GN model and 
there is a dynamical generation of a mass (gap) for the electrons. A global chemical potential /i is introduced in the 
theory to represent the extra electrons that are inserted in the system by the doping process. At a critical chemical 
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potential /ic = 0| [111, where Aq (« 0.7eV for irans-(CH)x) is the constant band-gap, the GN model undergoes a first 
order phase transition to a symmetry restored (zero gap) phase, agreeing with the experimentally observed non-metal 
to metal first-order transition in TPA [8,]. Although some times we shall refer to TPA throughout the paper (due to 
the large amount of data it has), the main results are applicable to any ID chain with the same structure of (CH)x. 

The critical doping concentration can be related with the critical chemical potential as yc = a/ie [3], where 
N = 2 can be regarded as an internal number (spin) of fermion flavors, h is the Planck's constant divided by 27r, 
vp = kph/ra is the Fermi velocity, kp is the Fermi wavenumber, m = j^toC? ^ and a (= 1.22A for irans-(CH)x)) is 
the lattice (equilibrium) spacing between the x coordinates of successive CH radicals in the undimerized structure, and 
io(« 3 eV for irans-(CH)x) is the intercarbon transfer integral for tt electrons, known as the "hopping parameter" . 
In Ref. [T^ the GN model has been used in the large- (mean-field) approximation as an effective model to describe 
the non-metal to metal phase transition in TPA, and they found a very good agreement with the experimentally 
measured yc(— 6%). 

In this work we use the GN model as a QFTM for doped (CII)x chains to investigate the consequences of a -Bo 
applied parallel to these "wires" at zero and finite temperature (T) . We show that the actuation of i?o produces an 
imbalance between the chemical potentials (and consequently in the densities) of the spin-up (f) and spin-down (|) 
conduction electrons inserted in the system. Depending on the level of imbalance between the chemical potentials of 
the two involved fermionic species the (CII)x chain can be in a non-condensed (zero gap) partially polarized state, or 
can be turned into a fully polarized state with a small but non-zero gap [T7| . It is shown that in irans-(CII)x with 
low imbalance (with respect to /ic), at zero T, the system acquires a partial spin polarization. With strong imbalance 
there is total absence of electrons from one of the two possible spin orientations (t or |) , meaning that the system is 
now fully polarized [t?] . We expect that this might have observable consequences in the magnetic properties of doped 
(CII)x or other conducting polymers. In the investigation of spin-polarized phases and the magnetic properties of ID 
cis-polymers, one has to consider the massive Gross- Neveu model [isj . as we shall do here. 

It is important to notice that the study of the effects of an external constant Zeeman magnetic field -Bq applied on 
a physical system is relevant in many physical situations, as in the investigation of metal-insulator transition [19| and 
magnetization [l^] in 2D electron systems. 

The paper is organized as follows. In Section II we introduce the discrete and continuum model Hamiltonians 
and the model Lagrangian describing (CH)x chains. In Section III the temperature dependent renormalized effective 
potential is obtained in the mean-field approximation and beyond. In this section we obtain an analytical expression 
for the effective potential at high temperature, as well as the chemical potentials dependent gap equation and the 
critical temperature at which the gap vanishes. Besides, the tricritical points of the second order transition line of the 
gap parameter are explicitly calculated by considering both homogeneous and inhomogeneous gap parameters. The 
temperature dependent magnetic properties of (CII)x chains are also obtained in this section. For completeness we 
study in Section IV the Massive Gross-Neveu model which is appropriate to describe Cis-Polyacetylene. In Section V 
we show a Summary of the main results and finally, we present the conclusions in Section VI. 

II. THE DISCRETE AND CONTINUUM MODELS 
A. The Model Hamiltonians 

1. The SSH Discrete Model Hamiltonian 
The tight-binding SSH Hamiltonian with electron-phonon interactions has the following form [lo| : 

H = + Hph + H^^^ph, (2.1) 

where 

= -ioXl('=n+i,sC„,s 4_,c„+i,s), (2.2) 
describes the tt— electron (e) hopping between site n and n+l. is the hopping constant for an undimerized structure, 

^p'^-E^ + fEK+i (2-3) 
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represents the lattice, i.e. the phonons (ph), Men is the mass of the CH group, K is the "spring" constant, u„ is the 
deviation from the undimerized structure on site n, and the interactions between these two is 

HT,-ph = a ^(u„+i - Un){cl+i^sCn,s + cl^sCn+l,s), (2.4) 

where a is the e — ph interaction constant. Summing and H-^^ph one defines in,n+i = U) — ct{un+i — Un)- 

2. Electron Correlations 

Electron correlations i.e., electron-electron interactions, are fundamental in the explanation of many physical phe- 
nomena displayed by conducting polymers as, for ex amp le, the observed negative spin densities in trans- (CH)x, which 
demonstrates electronic correlations in these chains [2l| . Since the SSH model (that is based only on electron-phonon 
interactions) gives zero spin density, other models that take into account electron correlations are naturally necessary. 
Besides, the inclusion of electron correlations are mandatory for the optical spectra of even Polyenes [23 |. 

The addition of a term describing electron interactions in the SSH Hamiltonian defines the Pariser-Parr-Pople 
(PPP) models [H, HI], which provides an extension of the Hubbard model by properly accounting for the long-range 
character of the electron-electron (Coulomb) repulsion. Thus, we have: 

Hppp = H + H-^^T^, (2.5) 
where H is the noninteracting Hamiltonian defined in Eq. (j2.ip and i?7r-7r is given by 

-ffTT-TT = + i^Vi,/'n;n/', (2.6) 

; 1,1' 

where U is an on-site Hubbard term, Vi^i' is an off-site interactions between electrons on site / and Z', and ni — 
^1 0*^1, a + cj pCi^p — 1 is the net charge density on site I with spin a {13). The prime in the last term implies I 7^ I' . The 
on-site correlations U, and inter-site Coulomb interactions Vi_i' are, in principle, arbitrary, and generally the Ohno's 
formula [Hi is employed [Ullil: 

Vu' = I ^ (2.7) 



1 



where ~ 1.29A, r; ;/ is the distance between sites I and V in units of A, and U — 11.13 — 11.26 eV. Thus, as pointed 
out in Ref. [23|, the geometry and U fix all inter-site interactions. 

The influence of electronic correlations on structural and electronic properties of irans-(CH)x, such as the ground- 
state energy, the amplitude of bond alternation and the effective force constant for the bond-stretching mode, has 
been investigated using the Hamiltonian of Eq. (|2.5p with Vi^v = but considering an effective U [2^. The study of 
Ref. [2^ finds a consistent modeling of trans-{CY{)x if U is chosen to be of the order of 7 — 9 eV. 

Since the consideration of electron correlations is out of the scope of the present work, in what follows we consider 
U^Vi^v =0. 

3. The TLM Continuum Model Hamiltonian 



As pointed out in Ref. [ll|, in the weak coupling limit, only electrons of the Fermi surface with energy of order w A 
(the gap energy) are affected by the dimerization. Then it is reasonable to linearize the band structure in the vicinity 
of the points ±/cf, where kp '^s, the Fermi momentum. After the linearization, TLM got the continuum version of 



H = — — ^- / dx 

2'KhvFXTLM 



A'{x) + 



A2 

"0 



dx ip'''^{x)[-ihvpa3dj; + A(a;)CTi]'0^(a;), (2.. 
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where the constants above are related with those of the original SSH Hamiltonian as hvp = "^toajXTLM = 
{2a)'^ /2TrKtQ, — AK/Mqyi is a dimensionless coupling, are the Pauli matrices, A(x) is a real gap related 

to lattice vibrations, and A(a;) = dA{x)/dt. ■0" is a two field component spinor — ( ] ^ representing the "left 



moving" and "right moving" electrons close to their Fermi energy, respectively, where s is an internal symmetry index 
(spin) that determines the effective degeneracy of the fermions, s = 1 =t, and s = 2 =1. 

B. The Model Lagrangian 

The Lagrangian density of the TLM model in the adiabatic approximation (A (a;) =0), is given by 

-Ctlm = V {ihdt - ihvFlA - 7oA(x)) r - TTT^ ^'(^) ' (2-9) 

s 

where the gamma matrices are given in terms of the Pauli matrices, as 70 = cti, and 75 = —(73. 
The four-fermion Lagrangian density of the massless GN model jl2| reads 

'Cgn = E - r + ^^1^ {i'^f , (2.10) 

where -0 — ip^jo and 'ipip = i/i^V* = X^fLi V'^V'^i ^-^d N is the number of flavors of Dirac fermions in the effective fleld 
theory. In the large-N (mean-fleld) approximation, one obtains [l^ 

£gn = V {thdt - iTiVFlsd.. - 7oA(a;)) r - A^(a:) . (2.11) 

Thus one sees that the equivalence between the TLM and the Gross-Neveu (GN) model is established by setting 
Atlm = iv^- Note that the fermion- lattice interaction in the equation above resembles the usual fermion-boson 



interaction in quantum field theory context 27 1 



In order to consider the application of an external Zeeman magnetic field Bq parallel to the system and its effects, 
it is convenient to start by writing the grand canonical partition function associated with /Cqn in the imaginary time 
formalism psl ]: 

Z^y" Di/i DV expjy dr J dx [icNlj , (2.12) 
where /3 = l/ksT, ks is the Boltzmann constant, and Lgn is the Euclidean GN Lagrangian density: 



Lgn = V il^'i-johdr + ihvFjA - A(a;) + 70/i + 705,5/15^0] t/-" - - — ^ A^{x), (2.13) 

where /x is the chemical potential, Ss — ±1/2, g is the effective g-factor and fiB — eh/2m « 5.788 x 10~^ eV T~^ 
is the Bohr magneton. The parallel magnetic field couples to the electrons' spin and produces the Zeeman splitting 
energy term AE = Ssg^J^sBo [H, |30| in Eq. (12.131) . From the form of the Zeeman energy term in Eq. (|2.13l) we see 
that it can be added to the chemical potential, thus defining an effective chemical potential terms in the Lagrangian 
density of the form, 

s=l,2 s=l,2 

= MtV'tVt + Mi^'iV; , (2-14) 

where = /i + S^, ^^ = 11 — and 5/i = ^iibBq- In [27j it has been chosen /i = /ic, such that when 5/i = Bq = 
the system is in the symmetry restored phase. Thus we rewrite Eq. (|2.13p as 
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■^GN = ^ Tp^i-johdr + ihvp^idx - A(a 

s=l,2 



GN 



(2.15) 



Integrating over the fermion fields leads to 



Z = exp l-^-^, / dx A^{x) \ n^^i detDj, 

[ nVFAGN J J 



(2.16) 



where Dj — — 7o9t + ifivp^idx + ToMj ^ ^i.'^) is the Dirac operator at finite temperature and density. Since A(a;) is 
static, one can transform Dj to the ujn plane, where uin — (2n + VjuT are the Matsubara frequencies for fermions, 
yielding Dj = (— ia;„ + /ij)7o + iTivrjidx — A.{x). After using an elementary identity ln{det{Dj)) — Trln(I?j), one 
can define the bare effective action for the static A (a;) condensate 



^e//[A] 
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r ^ 

/ dx A2(x) +^Trln(i:'j), 



(2.17) 



where the trace is to be taken over both Dirac and functional indices. The condition to find the stationary points of 
S'e//[A] reads 



5A{x) 



2/3 



hvpX' 



GN 



-A{x) 



6A{x) 



^Trln(I?,) 



(2.18) 



The equation above is a complicated and gene rally unknown functional equation for A(x), whose solution has been 
investigated at various times in the literature |3ll - l36| . Its solution is not only of academic interest, but has direct 
application in condensed matter physics as, for example, in (tI. [ssj. and in the present work. 

III. SYMMETRY BREAK AND CONDENSATES AT MEAN-FIELD AND BEYOND 

The GN model in ID has a discrete chiral symmetry which is dynamically broken by the non-perturbative vac- 
uum [31|]. As a result, there will be the generation of homogeneous and inhomogeneous A{x) condensates, as discussed 
below. 



A. Homogeneous A(a::) Condensates 



For a constant A field the Dirac operator reads Dj = (— ia;„ + /i,-)7n + ihvp^ip — A, so the trace in Eq. (j2.17p can 
be evaluated in a closed form for the asymmetrical ((5/i ^ 0) system [23|. From Eq. (|2.16p one obtains the "effective" 



potential Veff = InZ, where L is the length of the system: 



K//(A,Ait,;,T) 



fivwX 



-A^ 



F^GN 



dp 
2Trh 



2l3Ep + In ( 1 + e^'^^t ) + In ( 1 + e 



ln(l + e'^^")+ln(l + e-^^: 



where Ef^^^ = Ep± Ep = ^v\p^ + A^ 



-HE. 



(3.1) 



Zero Temperature and Zero Chemtcal Potentials 



The first term in the integration in p in Eq. (13. ip . corresponding to the vacuum part (/it,4. — T ~ 0), is divergent. 
Introducing a momentum cutoff A to regulate this part of Veff, one obtains, after renormalization where the mo- 
mentum cutoff is taken to infinity while maintaining finite and stable results, the following expression to the effective 
potential [ll,[i3|: 
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Tivp V-^GN ^TT/ TrnvF xntp J 

where mp is an arbitrary renormalization scale, with dimension of energy, introduced during the regularization process 
used to compute the appropriate momentum inte gral s. The minimization of Ve//(A) with respect to A gives Aq = 
and the well-known result for the non-trivial gap [12| : 

Ao = mpe^~ ^GN . (3.3) 

This is the phenomenon of dynamical mass generation in the massless model. From this equation it is easy to see that 
with the experimentally measured Aq and a, to and K which enters Atlm, one sets the value of mi^-. Equation (|3.2p 
can be expressed in a more convenient form in terms of Aq as 



VeffiA) 



ln(|)-l 



(3.4) 



which is clearly symmetric under A — >■ —A, that generates the discrete chiral symmetry of the GN model. As 
we pointed out before this discrete symmetry is dynamically broken and thus there is a kink solution interpolating 
between the two degenerate minima A = ±Ao of (|3.4I) at a; = ±oo [40| : 

A{x) = Aotanh(Aox). (3.5) 
In the next subsection we discuss the effects of space dependent A{x) condensates. 

2. Zero Temperature and Finite Chemical Potentials 
At finite chemical potentials and in the zero temperature limit, Eq. p.ip reads 

1,9 dp dp dp , , 



where p'^p} — ^y^/i^ ^ — A^ is the Fermi momentum of the spin-t(|) moving electron. We integrate in p, observing 
that the renormalization is the same as before, to obtain (T7| : 

Fe//(A,/it4) - K//(A) + eiJi + 92^2, (3.7) 
where T4//(A) is given by Eq. p.4p . Qi^2 = ©(m^^ ~ A^) is the step function, defined as 9(x) ~ 0, for a; < 0, and 
8(a;) = 1, for a; > 0, and 3'i,2 ^ ^ 



AMn( ^^'^^^?^ .l-Mt.,i./^;,-A^ 



Although finite chemical potential is being considered now, the effective potential is unaffected if both fi^ and fj.^ 
are smaller than A, in which case it is reduced to that of Eq. p. 41) . For /i-f-^j, > A the ground state is determined by 
jointly finding Ao(/i-|-,4,), and the analyzing of the effective potential at the minimum, Ve//(Ao(/x-|-4)). Extremizing 
Ve//(A, /i-|-_4,) with respect to A yields the trivial solution (A = 0), and 



where Si. 2 — In ^^^^ — ^a^'^ ^^"^ have made use of Eq. (13.31) to eliminate mp in the gap equation above. 

Equation (|3.8I) can not be solved in a closed form for A as a function of Nevertheless, for Bq = i.e., for (5/i = 
which implies /i-j- = /.tj, = fi, it can be easily solved: 



7 



Ao(m) = VAo(2a*- Ao) , (3.9) 

where Aq is given by Eq. p. 31) . and < /i < Aq. However, the analysis of Ve//(Ao(/i)) shows that Ao(/i) is a 
local maximum and the solution A — Aq represents a global minimum of 14// (A, /i) while /i < /^c, (where the critical 
chemical potential, fic = is obtained by Ve//(A = 0,/ic) = T4//(A = Aq, /ic)), with Aq — representing a false 
minimum. When fi > fj,c, A(u > jic) represents a local maximum and Aq = is turned into a minimum through a 
first order phase transition [isl, flGj. agreeing with experiment Q. Thus, Aq as a function of the chemical potential in 
the absence of Bo has the following expression: 

Ao(m) = e(/ic-A*)Ao. (3.10) 

The situation changes considerably when Bq is turned on, where A = still represents the minimum of 
Ve//(A, /i-|-.|) for < = 0.38Ao, where S^c is obtained from the equality Ve//(A ~ 0, /x-|-4((5^c)) = ''4//(A = 
Ao((5/ic), fJ-t-li^l^c)), with Ao((5/ic) being the solution of Eq. p.8p for ii^_i{dfic)- For ^/i > S/ic, the 02 function prevents 
the "4," term in the effective potential of participating in the minimum. In this case, Eq. p.8p can be rewritten as 

A^ - 2/itAgA + A;^ = 0. (3.11) 

3. Magnetic Properties at Zero Temperature 

As mentioned before, the application of a static magnetic field on the system results in a Zeeman energy given 
by /S.E = SzgiJ-BBo, where Sz = ±1/2, g (~ 2 for irans-(CH)x) is the effective 5— factor, fis = eT%/2m « 5.788 x 
10~^ eV T~^ is the Bohr magneton, m is the bare electron mass, and Bq is the magnetic field strength. Then we 
have 5^1 = \AE\ = ^^ibBq. 

The number densities n^^i = — ^ 14// (A, /^f.j.) are obviously imbalanced due to the asymmetry between /i-j- and 
fil, and will depend on Sfi. Before the critical asymmetry doping, Ao((5/i) — and the densities read 

n^,l{6fi/ Ao < 0.38) = ^^Mt,;- (3-12) 

For such a low imbalance, compared to the critical chemical potential asymmetry Sfic, the total number density is the 
same as in the symmetric limit or, in other words, is independent of the applied magnetic field: 

2 

riT = n+ + ri] = 2n = — /ic. (3.13) 

TTnVF 

In spite of that, there is a partial spin polarization 

2 2 

Sn = - ni ^ — — S^i = — — /is-Bo- (3-14) 
Trhvp Trnvp 

The zero T Pauli magnetization of the polarized chain is defined as: 

M = fiBSn= -P-Bq. (3.15) 

nnvp 

The magnetic susceptibility is given by 

oBq Trnvp 

which is the well known Pauli expression of the magnetic susceptibility for noninteracting electrons at zero temperature. 
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Increasing the asymmetry beyond the critical value Sfic = f /^sSo,c, we find 



nt(5M/Ao > 0.38) = y^M^- A2, (3.17) 
where A in the equation above is the solution of Eq. p. lip for a given S/jL > and 

7i^((5/i/Ao > 0.38) = 0, (3.18) 

meaning that, effectively, there are only spin-f electrons in the system i.e., the chain is fully polarized. The magneti- 
zation is now 



^/P = MBnt=^7;4^, (3.19) 

which also depends on the external magnetic field Bq through A. The results for the zero temperature magnetization 
of a (CH)x wire under i?o are displayed in Fig. [TJ 



Bp=o — — o > < o — ^-e — — <-e — e-> — *^ 



Bp<B(, — e+ — Q > o > — — e+ — e+ — <^ — 



Bq^Bq — — o > o > — — — — ^ — Q^^ — 

FIG. 1: Schematic representation of the spin orientations in the three possible situations in respect to the actuation of a Zeeman 
external constant magnetic field parallel to a (CH)x wire. In the first line _B = and consequently M = 0. In the second line 
B < Bo,c ~ 2SiJ.c/gf^B implying in a net magnetization of the wire. Finally, when B > _Bo,c the wire is fully magnetized. 



4- Finite Temperature and Chemical Potentials 

The effects of finite temperature and densit y h ave been taken into account in the symmetric GN model (5/i = 0) 
by several authors. See, for instance, refs. [l5l. l4ll. l42j and references in there. 
We can rewrite ^//(A, /^-|-,4., T) as 

Fe//(A,Ait,i>7^) - K//(A) +K//(m>7^), (3-20) 

where 

K,,(A.t,,,T) = -ksT J^^ S + + (l + ^""""O + + ^^"''0 + + ^'"""O 1 ■ (3-2^) 

Since we can not calculate expression (j3.2ip in a closed form, we shall use a high temperature expansion to evaluate 
it. Using the function 

/(a, 6) ^ dx [in (l + 6-^"'+"'-'') + In (l + e-^"'+"'+'')] , (3.22) 
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where a = A/ksT, and b = u/knT, which can be expanded in the high temperature hmit, a << 1 and b << 1, 
yielding, up to order a'' and 6^ [4l |. 



/(a«M«l) = - + ---ln(-)--(l 



7Bj 



7^) 
87r2 



i||^5V + 0(aV), (3.23) 



where je ~ 0.577... is the Euler constant and ^(n) is the Riemann zeta function, having the vahies ^(3) « 1.202, 
and ^(5) « 1.037. With the equation above, together with Eq. p.2p . the high temperature asymmetrical effective 
potential is written as 



A2 



In 



7b 



Ao 
7g(3) A^ 
87r2 (fcBT)2 



3hvF 



{kBTf 



(3.24) 



(fci3T)2 



64^4 (A^t + A'i) 



A^ 



The equation above may be rearranged in the form of a Ginzburg-Landau (GL) expansion of the grand potential 
density, which is appropriate to the analysis of the phase diagram in the region near the tricritical point. 



where 



ao(m'^) 
a4(Ait.i'^) 



2nhvf 



[f4 + f4] 



3hvp 



(3.25) 



(3.26) 



irfiVF 



In 



TTkBT\ 7^(3) (a^? + /^J) 



e'^-BAo 



i2T:^hvF{kBTf 



8^2 [kBTf 

186C(5) (Mt + m|) 



-14C(3) 



47r2 (fcBT)2 



(3.27) 



Extremizing Veff we find the trivial solution (A — 0) and the chemical potential and temperature dependent gap 
equation 



a2 



(3.28) 



which has meaning only if the ratio ^ is negative. Besides, a stable configuration (i.e., bounded from below) requires, 
up to this order, Q4 > 0. At the minimum Veff reads 



Veff,min — ClQ 



4^4' 



(3.29) 



The critical temperature Tc is the temperature at which the coefficient of A'^ vanishes i.e. a2 — 0, or 



In 



e^-Ao j 8^2 ^keTef 



0, 



(3.30) 



where Tc — Tc{fi-\-, /ij.). As will become clear below, the equation above defines a second-order transition line separating 
the non-metallic (A 7^ 0) and metallic phases (A = 0). At /i^ = /i^ = 0, we recover the well-known result for the 
temperature at which the discrete chiral symmetry is restored [46i |: 



TM = ^4 = 0) = Te{0) 



Ao 

TT kB 



(3.31) 
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In order to find Tc{fi-\- , fii) we define dimensionless variables v = "^^^^i (1^7^(0))^ ^ ~ t'\o) ' ^^'^ with tlie help of 
Eq. ([33T|) we rewrite the L.H.S. of Eq. ([330l) as 

y{t) = ln(i) + ^. (3.32) 

The zeros of y{t) for a given ly, i.e., for a given /i^ + /ij, are the respective T^. This defines the (second-order) 
versus /x| + /ij phase diagram. A graphical inspection of y{t) shows that there is no solution for this function for 
v above certain value, that we define Vt^- Besides, at Vtc we have y ^ y' = 0. These two equations give ttc and z/f^ for 
the tricritical point Ptc = (vtcUc)- Strictly speaking, y and y' are associated with the coefficients of the second-order 
and forth-order terms of the effective potential expanded in powers of A [48] . Solving the equations y — and y' = 
self-consistently (which is equivalent to solve a2 = 04 = 0), we obtain the tricritical point analytically: 

_ 7^(3) (^t + Mj)tc _ 1 _ Ttc _ _ 1 

~ 8772 {ksTMy ~ 2^' 7^' ^ ' 

For v above certain value and less than vtc, the function y presents two solutions for Tc- However, the lower of these 
always corresponds to unstable solutions. The second-order transition curve, defined as the line starting at the point 
(0,Tc(0)) and ending at the point ((/i| + ^'^)tc,Ttc), comes simply from the solution of the gap equation. This curve, 
shown in Fig. [31 represents a system at finite temperature where the chemical potentials fi^ + = 2(/i^ + 611^) start 

from zero and increases until (/i^ + fJ^1)tc- Note that p,^ + Sfj,^ is zero if and only if p,^ and 6^^ are both zero. It is 
well known that below the tricritical point one has to properly minimize the effective potential rather than using the 
gap equation as the transition becomes first order. Thus Eq. (|3.30p cannot be used for finding Tc below Ptc since this 
equation is valid only for the second-order transition. In this case Tc has to be find numerically, through the equality 
Veff{Tc, ^^,ni,A = Amin) = Vcff{Tc, , fj.^, A ^ 0) , whcre A„iin is the non-trivial minimum of K//. 




FIG. 2: The phase diagram t — 
ends at the tricritical point Ptc - 



'^811-^^ (fcgT {0))'^ ' i^om Eq. (|3.30p . Tlie second order transition line 



jtjqj as a function of = 

- (^5 ■^)- Below this point the transition is of first order. 



The number densities 



^K//(r,Mt,i,A) read 



dp 



(3.34) 
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where nk{E^"^ ) 



e t,4- +1 



is the Fermi distribution function. The density difference 



6n 



(3.35) 



is zero if S^i = ^^ibBq — 0, at any temperature, since in this case we have the equahties nk{E^) = nk{E^), and 
nk{E^) — nk{E^). The physical meaning of these results is that at zero external Zeeman magnetic field, the f (up) 
and 4, (down) electrons of the conduction (+) band have the same density, and the same for the electrons of the 
valence (— ) band. 

In the high temperature regime, the number densities are given by 



irfivp 



1 



7^(3) A2 



47r2 {kBTf 



(3.36) 



In the high temperature limit the total number density, nxiT) — n^{T) + ni{T), is independent of the applied field, 
as happens at zero T 



nriT) = 



nhvp 



1 - 



7^(3) A2 



47r2 {kBTf 



(3.37) 



and for the density difference we obtain 



irfivp 



1 



7C(3) A2 



47r2 {kBTf 



(3.38) 



that, as we have observed before, is clearly zero if Bq — Sji ~ 0. Since the densities have to be evaluated at the 
minimum of the effective potential^ , we use Eq. p.28p in the equation above and find the temperatures at which the 
densities and, consequently, the density difference vanish. These temperatures are the solutions of 



7s 



In 



Ao 



TrkBT* 



7g(3) (m^ + 

8^2 (fcBT*)2 



0, 



(3.39) 



where T* — T*(/i-|-, /i|). As for the gap parameter, the equation above defines the second-order line for the densities 
and the density imbalance. At fi-f ~ fi^ — 0, we get 



7^c(A^t = Mi = 0)^T;(0) = 



Ao 



(3.40) 



It is very easy to see that T*(0) = which coincides with Ttc, where Ttc is given by Eq. p.33p . Proceeding as 

before we find 



7e(3)(A^? + A^!)*c _ 1 _ _ (3.41) 



87r2 (fcBr*(0))2 2e2' T*{0) 

defining the tricritical point for the densities, total density and density imbalance second order curves. 

Let us now verify the possibility of a fully polarized state at finite temperature. It would be possible with a magnetic 
field with a intensity such that in Eq. (13.36^ vanishes. In this case fi^ = fl — Sfj,c = fJ-c — f Ms^o.c = 0, or 

Bo.c-^. (3.42) 

5MB 



^ Remembering that the thermodynamical potential per volume is defined as the free energy density (or effective potential in the field 
theory language we are using here) at its minimum, Q{T, fi^^^) = V^f f{T, fif ^^, Amin)- 
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5. Magnetic Properties at High Temperature 
The Pauli magnetization of the chain in the high temperature hmit has the following expression: 



Mhigh t{T) = IJ-B^Tihigh t{T) 



In 



2MB 
■kTivp 

T 



1 



7^(3) 



47r2 [kBTf 
7^(3) (4 + 



87r2 {kBTf 



Bq, 



(3.43) 



The second-order 



where Tc(0) is given by Eq. (I3.3ip and we have made use of Eq. p. 281) to leading order in (j^^j'y 

line where the magnetization vanishes is the same as the one given by Eq. p.39p . With Mhigh t{T) we are in condition 
to obtain the magnetic susceptibility in this regime 



frr\ 9Mhigh t(T) 

XMgh t{T) = = x(0) - 



where 



X(0) 



Trhvp 



is the zero temperature Pauli magnetic susceptibility for noninteracting electrons given by Eq. p.l6p . and 



(3.44) 



(3.45) 



xiT) 



nhvF 



In 



T 



7^(3) 



T,(0); 47r2(fcBT)2 



(3.46) 



The function Xhigh t (T) also behaves at finite temperature as the densities and the magnetization, with a second-order 
transition up to a tricritical point given by Eq. p.4ip . Below this point the transition is again of first order. 

Note that despite the fully polarization, at Bq c the magnetization is given by exactly the same expression shown 
in Eq. ([Xig)) . With the help of Eq. ([XM]) we find: 



Mhigh T,c{T) = t^Bn^{T) 



jJ'B 

ttTivf 



1 - 



7^(3) A2 



47r2 (fcBr)2 
T 



In 



r*(o) 



{flc + SfJ-c) 

7^(3) (m^ + m!) " 
87r2 (fcsT)2 



0,c- 



(3.47) 



This shows that this function is indeed continuous for < Bq < Bq , 



B. Inhomogeneous A(a;) Condensates 

Since we consider the addition of a chemical potential (i.e., doping) in the theory representing a Peierls distorted 
material like TPA, and the effects of a Zeeman magnetic field on these materials, some important remarks are in order. 
It is well-known that doping in conducting polymers with degenerate ground states results in lattice deformation, or 
non-linear excitations, such as kink solitons and polarons, meaning that A(a;) can vary in space 0, H^iHJ- Therefore, 
one may expect not only homogeneous-like configurations (as considered in the previous subsection), but also that 
the inclusion of these excitations in any theoretical calculation in this model should be considered. In this context, 
within the GN field theory model that we are considering, by taking into account kink-like configurations in the large 
N approximation, the authors of Refs. [33| - (36| found evidence for a crystalline phase that shows up in the extreme 
T ~ and large fj, part of the phase diagram, while the other extreme of the phase diagram, for large T and small fi, 
seemed to remain identical to the usual large N results for the critical temperature and tricritical points, which are 
well-known results [lH| for the GN model. The restoration of the broken discrete chiral symmetry of the massless GN 
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model, which is signalized when the kinkantikink crystalline condensate [l^, |37 



is explained when an exotic supersymmetric structure is properly considered 39 1 



transforms into the kink crystal, 



To take into account the effects of inhomogeneous configurations in the GL expansion of the grand potential density, 
let us write Eq. p.25p in terms of A(a;) and its derivatives up to a4 [33l-[36|: 



(3.. 



where A(a;)' = dA{x)/dx. A straightforward variational calculation gives the following condition for the minimization 
of the free energy E = J Veff{x) dx: 



A{x)" - 2A{xf - —A{x) = 0. 



(3.49) 



The general solution of an equation of the form 



A{x)" - 2A{xf + (1 + ^)AlA{x) = 0, (3.50) 

can be written as [3a| 



A(a;,<j) = A(a:;) AoV^ sn(Aoa;;<j), (3.51) 

where sn is the Jacobi elliptic function with the real elliptic parameter < <7 < 1. The sn function has period 
2K(<;), where K(<;) = Jq^^[1 — sin^(z)]~-^/^dz is the complete elliptic integral of first kind. At T = and finite 
density (or /i), the <; parameter can be determined by a variational procedure yielding an average fermion density 
dependent transcendental equation for which can be solved in the low- and high-density limits [35]. At finite T and 
/i, the <; parameter can be determined within the Schnetz, Thies and Urlichs' ansatz to solve the inhomogeneous gap 
parameter problem, together with the scale factor ^ by a numerical minimization of the renormalized grand canonical 
potential [33]. 

A(a;) in (|3.5ip represents an array of real kinks. When = 1 Eq. p.5ip is reduced to the single kink condensate 
given in Eq. p.5p . By a direct comparison between Eqs. p.49p and (I3.50p one can easily identify the scale parameter 
Ao at high T (i.e. T near as 



A„^^(-^)^.A„(.). (3.52) 

Since > 0, the solution for inhomogeneous condensates has physical meaning only if the ratio ^ is negative, as 
in the case of homogeneous condensates. 

In order to see the effects of Eq. (I3.51|) in the grand potential density, we add the term —a4^{A{x)'^)" to Veff{x) in 
Eq. p.48p . It is easy to see that this term does not alter the solution for the minimum of the free energy, Eq. p.49p . 
but allows the use of the following important identities supported by the solution in p.50p : 



(AixYf = A\x)-il+<;)AlA{xf+^A*, (3.53) 
iA{xf)" = 6A4(x) -4(l + 0A2A(a;)' +2^A^. 

With the equations above, we can write the x-dependent grand potential density as 



Veffix) =ao + a2Aixy + a^- [(1 + ^)AlA{xf + ^A^] . (3.54) 

Averaging over one period, it is found [s^ that < A(x)^ ^1 — Aq, where E(<;) is the complete elliptic 

integral of second kind. The ratio E((^)/K(<r) is a smooth function of interpolating monotonically between and 1. 
Thus we can write 



14 



<Veff[x)>=ao+A2/:^l+Ai/S.l, (3.55) 

where 

A2 = a2 f 1 - ^) , (3.56) 



A4 = Q!4- 



K(0 

^ + (1 + 



E(0 
K(0 



As we will see next, with the consideration of inhomogeneous A(x) condensates we obtain very interesting results. 
This generalizes the analysis of the grand potential density, such that the results obtained in subsection A for 
homogeneous condensates are particular cases which are obtained when specific values of the elliptic parameter ^ are 
taken: 

(I.) For = 0, ^1^"°! = 1, which implies A2 = A4 = 0, and the grand potential density is that of the 
metallic phase, for which A = and 14// = Q;o- 

(II.) For = 1, = 0, so the grand potential density is that of the homogeneous case, Eq. (|3.25p . at 

the non-trivial minimum Aq: Veff = ao + (^2^0 + a4^Q- 

Besides, we find that: 

(III.) For a < <r < 1, the grand potential density is that of Eq. p.55p that, as expected, is the interpola- 
tion of the two previous ones. 

And finally, we find that the tricritical point is "stationary" under the application of an external Zeeman 
magnetic field on CH chains: 

(IV.) For any < ? < 1 it is easy to see from Eq. p.55p that the tricritical point is still found for a2 = ct4 — 0. 
These coefficients are fi-\-^i and T dependent and were not affected by the space dependence of the condensate A{x), 
as seen in the symmetric case (i5/i = = 0) [33l - [36j . The novelty is that the location of the tricritical point of the 
TPA phase diagram under the influence of an external Zeeman magnetic field is also unaltered even considering a x 
dependent grand potential density. This happens because in the high temperature limit the influence of the Zeeman 
field is not sufficient to change the position of the tricritical point. 

C. The Gross-Neveu Model Beyond the Mean-Field Approximation 

We now investigate the GN model beyond the mean field or large N approximation employing the optimized 
perturbation theory (OPT). The OPT is a method to take into account non-perturbative contributions to the mean- 
field results, which are of crucial importance to recover the reliability of perturbation theory at finite temperature, 
mainly in theories displaying thermodynamical phase transitions (43j . For a better description of the implementation 
of the OPT procedure, please see refs. [4l,|4§|. The OPT is closely related with the optimized expansion (OE) method. 
As OPT, the OE is an alternative non-perturbative technique which also allows one to go beyond the framework of 
the standard perturbative expansion [45]. 

Applying the usual OPT interpolation prescription to the original Lagrangian density given by Eq. (|2.10p . we define 
the interpolated theory: 

^gM,4') = ^V^' i^^lodt - invF-fA)r - ?7(1 - + S^hvFii'^f ■ (3.57) 



It is easy to verify that at S ^ we have a theory of free fermions, and the original theory is recovered for S — 1. 
•q is an arbitrary mass parameter. Now, by re-introducing the scalar field A, which can be achieved by adding the 
quadratic term corresponding to a Hubbard-Stratonovich trick: 
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5N 



GN 



(3.58) 



to Lqt^[iIj^ ijj), one obtains the interpolated model corresponding to the original TLM model given by Eq. (|2.9p 



s=l,2 



8N 



2hvp\Gi 



(3.59) 



Introducing the chemical potentials as before and considering the actuation of the Zeeman field, as explained below 
Eq. ((2?T2t . we write: 



5N 



s=l,2 



2hvpX 



GN 



(3.60) 



Since Eq. p.60p in the limit (5/^ — > is the same model already studied in Ref. [4J| , so we do not repeat all the details 
related to the free energy density derivation here, where only the main steps and results relevant for our application 
to the CH-chains will be presented. 

Generally, the OPT method can be implemented as follows. Any physical quantity, F'^'^\ is perturbatively computed 
from the interpolated model, up to some finite order-fc in S, which is formally used only as a bookkeeping parameter 
which is set to the unity, at the end of calculation. But in this process any (perturbative) result at order k in the OPT 
remains rj dependent. This arbitrary (a priori) parameter is then fixed by a variational method that then generates 
nonperturbative results, in the sense that it resums to all orders a certain class of perturbative contributions through 
self-consistent equations. Such optimization method is known as the principle of minimal sensitivity (PMS) and 
amounts to require that F^''^ be evaluated at the point where it is less sensitive to this parameter. This criterion 
translates into the variational relation ISOl 



dr] 



5—l.ri—rj 



= 



(3.61) 



The optimum value fj that satisfies Eq. (j3.61l) must be a function of the original parameters, including the couplings, 
thus generating "non-perturbative" results. In our case, we are interested in evaluating the optimized free energy at 
finite temperature and density for the scalar field, A, once the fermions have been integrated out. 



1. The Optimized Free Energy Density 

We generalize to the asymmetrical case the Landau's free energy density (or effective potential, in the language 
of quantum field theories) to order-5, that has been evaluated in Ref. j43 | using usual functional and diagrammatic 
techniques. The result is: 



J(A,7?,T,Ai^,;) = 5 



iVA2 



N 



2Xq^vfTT' 2'KVp'h 



rj 



2 \ ri 



+ (fcT)2[/i(,;,^^,r)+/i(,7,/.;,T)] 



Mviv-A) 



TTVph 
g ^GN 2 



In ( (12(^7, Mt,r) + /2(r/,Mi,r)) 



Inf^ 



1 



i(fcT)2[/f(r;,^,,T) + /|(^,^,,r)] 



An^vph 2 

where the functions /i , I2 and I3 are given respectively by 



(3.62) 



Ii{ri,fi,T) = J dx |ln 



1 , 1 \ ( ^ '\ M 



In 



l + e-v/^^^(¥f+if^' 



(3.63) 
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dx 



y^^TWIW LeV^^^(^+w + 1 eV-^+(*)'-* + 1. 



(3.64) 



and 



sinh I 

V/cT 



cosh 



(fcTp 



^ +cosh(^) 



(3.65) 



In Eq. p.62p . A is the homogeneous (constant field) configuration for the scalar field discussed already in the sub- 
section IIII Al In the computation performed in Ref. 0| , the free energy density has been renormalized using the 
MS scheme for dimensional regularization. To a direct comparison between our results and that of Ref. 4J] , we have 
adopted dimensional regularization in this section instead of the cutoff regularization we employed in the previous 
section. As will be clear soon, Eq. p.62p . evaluated at first order in the OPT, already takes into account corrections 
beyond the large- (or mean-field) leading order result. 

By optimizing Eq. p. 621) through the PMS condition, Eq. (I3.6ip . wc obtain the optimum value for the parameter 
1] ~ f], which is then re- inserted back in Eq. (j3.62p . allowing us to compute the order parameter A that minimizes 
the free energy. Using the PMS procedure we then obtain, from Eq. p.62p a.t rj = fj and (5=1, the general result that 
factors into: 



^7- A + 77^^y(?7,/xt4,T) 



AttN 

where we have defined the function 



= 0, 



(3.66) 



y(r?,Mt,i:r) = ln( -i(/2(r;, A.^,T) + 12(77, A^i,T)) 



(3.67) 



Note that when A^ — > 00 in Eg. (13.661) . one obtains the solution ij — A and the mean- field standard result is exactly 
reproduced, as it should be (44ll5l|. 



2. The Gap Energy at Zero Temperature and Chemical Potential Beyond Mean-Field 

In order to perform a numerical analysis we must fix all parameters. This can be done by considering the gap 
energy. In the GN language the order parameter A is just the TLM gap parameter which, at T = and /J. = 0, has 
been denoted as Ao in Eq. At T = ^ = we have that h = I2 = h = so that Eq. ([Slee]) and dJ/dA = 

give, at order-(5, respectively: 



A = 77 



Agn / wf 
27rA^ ^\f] 



(3.68) 



and 



. Agn / mp 
A = 77 in 

TT V 77 



(3.69) 



Plugging 77 In (^^^ = ^ from Eq. (ISTM)) in Eq. ([5:^ we find. 



A= 1- — 
2A^ 



(3.70) 
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Solving these equations self-consistently yields: 

f]si (T = /^t,i = 0) = %i (0) = mpe ^ (3,71) 

and 



Asi{T = = 0) = Asi{0) = F{N)-^mFe ^g/*-(«) ^ 



(3.72) 



where F{N) = (l — 2^). As pointed out in Ref. Q, Eq. (I3.72p exphcitly includes corrections beyond large N, 
as obtained from the OPT approach, to the mean-field result. More precisely, taking the mean-field approximation, 
TV — 00 in Eq. (|3.72p and using the relation Aqn = ^""Atlm, the OPT result exactly recovers the mean field result 
for = 2 as expected 0, [5l|. Besides, we have that A5i(0) in the limit N 00 reproduces Eq. p.72p with an 
appropriate redefinition of the (arbitrary) renormalization scale mp- 

It is worth to mention that, usually, in a renormalizable quantum field theory, one can choose arbitrary value for 
nip and Agn will run with the scale appropriately, at a given perturbative order, so that Aq remains scale-invariant 
as dictated by the renormalization group ^ . 

3. The Gap Energy at Zero Temperature and Finite Chemical Potential Beyond Mean-Field 



We start by taking the T -H- limit of the functions defined by Ens. ((3^ . ((3^ and l^M\i which appear in the 
free energy density, Eq. p.62p . In this limit, these equations read [44|: 



\iuiT'hir^,^l,F)^--e{p~v) 

T~¥0 I 



if In 



lim /2(?7, ^JL, F) = 9{^jL - -q) In 



lim F^l3{'q,fi,F) = sgn(/i)6'(/i - r/) _ ^2 



(3.73) 

(3.74) 
(3.75) 



Using Eqs. (I3.73p . p.74p and (|3.75l) in the PMS equation, Eq. (|3.66l) . we notice that it can be divided into two 
cases: i) /if^i < rj, which is equivalent to the situation where F = /i^,^ = studied in the previous subsection, and ii) 
Mt,-I- ^ V- For ^f^i > rj we have 



27rA^ 



where 



27rA 



(3.76) 



y (77, T = 0) ^ y (r;, /.t,;) = In (^^) - ^ 



In 



In 



(3.77) 



and we have defined the function 



G(77,m) = G 



Mr 



(3.78) 



Using 



(3.79) 
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from dS^/dA = in Eq. (|3.76p . we find a second order equation for A: 



2N 



'^GN 



whose solution is: 



(3.80) 



A = - 



2(1-2^) 



i.^ri-^iG 



1 

27V 



(3.81) 



We have taken the positive solution from Eq. p.80p since the limit A — )■ of Eq. p.8ip correctly reproduces the first 
order in S solution displayed in Eq. p.70p . The equation above can be easily expanded in powers of the coupling 
constant A^jv: 



(l-w)l 27V V 2iV 



GN 



2N \ 2N \ TT \ TT 



(3.82) 



We can plug A from Eq. (I3.79P in Eq. (I3.82p and obtain an equation for ry which, to small values of , reads 

TT 1 



y('7,/^t 



I) 



(3.83) 



Xgn (1 - 27v) ' 

where ^ is given by Eq. (j3.77p . The equation above can not be solved in a closed form for fj as a function of A^t.J-- 
However, for Sfi — i.e., for fif = /i| = /i it can be easily solved and generalizes Eq. p. Tip for non-zero values of the 
chemical potential, yielding: 



%i(m) = V%i(0)(2Ai-%i(0)) 



where 7751(0) is given by Eq. p.7ip . and ^^^^^ < H < ?y5i(0). With this equation we find to this order 

1 



Asiifi) = 1- 



2N 



(3.84) 



(3.85) 



As for Agi (0) given by Eq. p.72p . Eq. (|3.85p obtained from the OPT approach, explicitly includes corrections beyond 
the large N to the mean-field result. Note that the equation above is reduced to the mean- field (local maximum) 
result, given by Eq. (|3.9p . in the limit N ^ 00, as it should. 



4- The Number Densities at Zero Temperature Beyond Mean-Fteld 

As we have seen before, the densities are obtained by the usual relation 71-1-4 = — ^ Vef / {T, /i-f.^. A) at A = 

Amin = A. Now since we are considering corrections to the large- results, we have to include the other minimum 
condition = fj. Then we want to know how the corrections beyond the mean field will affect the densities given by 
Eq. p.34p . Then one obtains 



nt4(T,Mt,Mi) = ^ 



N 



Ztt 



A)^I^{v,f^t,i:T) 
Ztt 



•agn 
4^ 



(3.86) 

jj=,-(.A=A 

where the primes indicate derivatives with respect to /if, 4.. 

Since in this subsection we are interested in obtain the zero temperature corrections to the densities beyond mean- 
field, the thermal effects will be neglected. 

Using Eqs. ((3J3)) . ((3J4)) and (|3J5)) in Landau's free energy density, Eq. (|3.62p . we observe that it can again 
be divided into two cases: i) /x^^ < ry and ii) ^-1-4 > 77. The case /z-f- j^ < 77 is equivalent to the situation where 
T — 4, — giving 7t,-|-.4_(T — Q,^^ < 7/,/i^ < 7/) = 0. For case ii) one has to resort to the numerical solution of the 
nonlinear equations above. 
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IV. THE MASSIVE GROSS-NEVEU MODEL AND THE CIS-POLYACETYLENE 

A. The MGN Model Lagrangian 

The cis-polyacetylene has non-degenerate ground states and is appropriately described by the massive Gross-Neveu 
model (MGN) model. The MGN model is obtained when the discrete chiral symmetry is explicitly broken by adding 
to the Lagrangian (|2.13p a fermion bare mass term M for both the spin-^ and spin-| electrons: 



3 = 1,2 



+ ihvFjid.j; - A{x) - M + 7o/ij]i/'j 



1 



(4.1) 



The bare mass term M corresponds [3J] to the Brazovskii and Kirova's external contribution from the rigid polymer 
skeleton to the total gap [ij]. As we will see below, the presence of a (bare) mass term ensures the elimination of 
kink-antikink configurations, which are suppressed in the thermodynamical limit even in presence of finite chemical 
potentials for the electrons 54]. Then in this section we examine only homogeneous condensates. In presence of a 
mass term it is more convenient to redefine the A field by shifting it by a constant A — A + M. Then the Lagrangian 
(apart from constant terms) reads: 



MGN 



E 

i=i,2 



ipj[--iofidr + ihvFlidx - A + 7oMj]'i/'i 



2M 



-A- 



1 



-A^ 



(4.2) 



1. The Effective Potential of the MGN Model at Zero Temperature and Zero Chemical Potentials 

Proceeding as before, we obtain the following expression for the renormalized effective potential of the MGN model 
at zero temperature and chemical potentials 



A2 / 1 



.^InfA 



Substituting mp from the solution for A when M = we find: 



2M 



-A. 



(4.3) 



K//(A,M) = 
The minimization for Vef f{A., M) gives: 



2n'hvF 



-4 



- 1 



2M 



-A. 



(4.4) 



/A2 



27rM 



0. 



(4.5) 



^0/ ^GN 

From the equation above it is easy to see that when M = the solutions are A = and A ~ ±Ao at a; = ±oo, 

interpolated by a kink, as seen in Section [nil Besides, we have that Ve//(A = 0) = and Ve//(A = ±Ao) = — ■ 

Then, since A = corresponds to a local maximum, the minima is degenerate at ±mFe^ ^"^n . For any M ^ the 
degenerate vacuum (minima) disappear, as well as the kink, giving place to a unique solution for a given value of 
M [5^. As pointed out in Ref. (5^ if any soliton exists at all, its stability has to depend on the energetics of trapping 
fermions. 

Eq. (I4.5P is very sensitive to the value of M. In other words, there will be solutions for (positive) A only up to a 
maximum value of M. If we define a dimensionless function f{x, M) = f{x) ~ x^ In (a;^) + Mx, where cc = ^ and 



M = 1—^ 
Ao(M): 



we have to solve self-consistently f{x) — and /'(a 



^ 4f(£) ^ 
dx 



to obtain the windows for M and 



0< M < 



AgnAq 



— < Ao(M) < Ao. 
e 



The effect of a mass term in the MGN model can be summarized simply as: as M is increased from to M„ 
the minimum of 14// (A, M) moves smoothly from Aq to 



(4.6) 



7re ' 
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2. The Elective Potential of the MGN Model at Zero Temperature and Finite Chemical Potentials 

Since after the change of variable A — > A + A/ the mass parameter M enters only in the classical part of the effective 
potential, it can be written as 

ye//(A,M,Mt4) = K//(A,Af) + eiJi +62^2, (4.7) 

where Vef f{^,M) is given by Eq. (|4.4p and Ji^2 are defined after Eq. (|3.7p . The extremization of T4//(A, M, 
yields 



27rA/ 
Agn 



0, (4.8) 



where gi^2 = hi ( ^'^'^'''^^'^ -^ ) . In the hmit M Eq. reduces to Eq. ((5^ . The equation above can be 



solved only numerically for arbitrary values of /^^^^ and AI . However, for (5/i = we obtain the following fj, and M 
dependent equation to first order in AI/A: 

A3 + (A2-2AoAi)A-^(Ao-A^)Ao=0. (4.9) 
Agn 



This equation has the analytical solutions: 

1. When/i^AoA A = Ao(^)'^'; 

2. When /i — Aq, A = 0. The other solution, A — Aq, is a local maximum. This can be understood simply as: for 
small AI and such a large value of /x the broken symmetry has been restored. 

For other values of /i, Eq. (I4.9p has to be solved numerically. 



V. SUMMARY OF THE RESULTS 



We show in Tab. |T] several physical quantities calculated within the GN model at mean-field and OPT (beyond 
mean- field) approximations at various medium conditions, i. e., temperature and densities. The places in the table 
with the ? means that, to the best of our knowledge, there are no studies concerning that physical quantities at 
mean-field and beyond in that situations. 



VI. CONCLUSIONS 



We have studied the influence of an external Zeeman magnetic field Bq on ID (CH)x chains, which can be converted 
into a partially or fully polarized organic conductor, depending on the intensity of the Bq applied on the wire. The 
mean-field finite temperature phase diagram of the chain has been investigated. We found that the gap parameter 
and the magnetization (as well as the densities and density imbalance, and the magnetic susceptibility) of this ID 
material under Bq have a similar second-order behavior until their respective tricritical points are reached. Below 
these points the transitions are of first order. We found these two tricritical points analytically. 

He have shown that the location of the tricritical point in the t yJo}^ versus v '^^^^^ (k^'r'io))^ ) phase 
diagram stays at the same place by considering both homogeneous and inhomogeneous condensates, as occur in 
symmetric = Bq = systems |33l - [3a |. We have found the critical magnetic field Bq c necessary to have a fully 
polarized organic conductor at zero and finite temperature. We found that for a given magnetic field below i?o.c, 
partial polarizations (magnetizations) can be realized experimentally, provided the temperatures are kept outside the 
"non — metallid'' region of Fig. [51 It is worth noting that, according to Eq. (|3.42p . for other ID systems with a smaller 
critical chemical potential or with a greater effective (^-factor, a smaller critical magnetic field necessary for a fully 
polarization of these ID chains would be found. 

In Table U we show the homogeneous and inhomogeneous gap parameter A{x), the Pauli magnetization AI and 
the magnetic susceptibility x, and zero and finite T and /i calculated at mean-field approximation and beyond. The 
quantities identified with the ? mean that further investigation are needed to find them. 
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Approximation \ Physical quantity 


Homogeneous A 


Inhomogeneous A(2;) 


M 


X 


MF at zero fi and zero T 
MF at finite /i and zero T 


Z\q TTLpC -^GN 


* 


u 



u 



IVIF at finite /i-]- and zero T 


Num. Solut. of Eq. ((5^ 


? 


Bo 

irhv f ' 




A/TP finifo 11^ , QTirl ViinVi T~' 


/ Q12 
V 2a4 




-CjQ. l|o.4o|) 




Beyond MF at zero T and zero fi 


F(iV)~imFe ^GN^-(") 


? 


? 


? 


Beyond MF at finite fi-^^^ and zero T 


Num. Solut. of Eq. (|3.83|) 


? 


? 


? 


Beyond MF at finite fi^^^ and T 


? 


? 


? 


? 



TABLE I: Physical quantities calculated at MF and beyond approximations and situations for the temperature and chemical 
potentials. The inhomogeneous gap parameter at zero T and finite /i marked with the * symbol is discussed below Eq. 13.511 
Ao(<;) is given by Eq. (|3.52|) . F{N) = 1 — 1/2N, M is the Pauli magnetization and x is the magnetic susceptibility which are, 
of course, zero for (5^ = |^sBo = 0. 



We have also seen that the consideration of electron correlations in the investigation of some specific properties of 
ID CH chains is of crucial importance. Then we intent to analyze a continuum version of the Hamiltonian of Eq. (j2.6p 
within the quantum field theory approach, and publish it elsewhere. 

A direct extension of the present work would be the investigation of the possibility of using polymers or graphene 
ribons as spin-polarized conductors, upon the joint application of a Zeeman and an electric field in the system. This 
would enable to study the transport properties of ID chains at zero and finite temperature. 

As a next step in the investigation of the effects of magnetic fields on ID systems, we would like to study the results 
of parallel and perpendicular magnetic fields applied simultaneously and independently on a polymer or graphene 
wire. This could open interesting possibilities for using these type of materials in practical applications as electronic 
devices. 

The methods we used here, based on the calculation of the (grand) effective potential free energy, which are more 
suitable for the analysis of the phase structure of the model, are useful also for the study of the consequences of a 
perpendicular magnetic field applied in the system. Indeed, one can incorporate the effects of a magnetic background 
on the effective potential and investigate how the chiral transition and the condensate (in our case the A gap) are 
modified. Assuming that the system is now in the presence of a strong magnetic field background that is constant and 
homogeneous, one can compute the modified effective potential following the procedure presented in detail in Refs. 
[5a459| . We hope to make further studies on the magnetic properties of these systems based on the model we have 
studied here, and to report on them in the future. 
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